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f~* > ; ■ We investigate dynamics of the homogeneous time-dependent SU(2) 

' Yang-Mills fields governed by the non-Abelian Born-Infeld lagrangian 

which arises in superstring theory as a result of summation of all or- 
^5 ders in the string slope parameter a'. It is shown that generically the 

Born-Infeld dynamics is less chaotic than that in the ordinary Yang-Mills 
theory, and at high enough field strength the Yang-Mills chaos is stabi- 
^ lized. More generally, a smothering effect of the string non-locality on 

(-H , behavior of classical fields is conjectured. 

The chaotic nature of the classical Yang-Mills dynamics was recognized more 
than twenty years ago ^[2IEIEIE1- It can be revealed already on simple models 
of the homogeneous fields depending only on time Further studies have 
shown persistence of chaotic features in classical dynamics of spatially varying 
Yang-Mills fields leading to a new kind of the "ultraviolet catastrophe" [7|. 

Here we would like to investigate whether this tendency is preserved on a 
deeper level of superstring theory, or the latter provides a smothering effect. 
The possibility to explore this question quantitatively is related to existence of 
a closed form effective action for the open strings or the stuck of D-branes which 
accumulates all orders in the string slope parameter a'. Strictly speaking, an 
exact in a' effective action is known only for the U(l) gauge field, in which 
case it is the Born-Infeld (BI) action |SJ |U] . In the non-Abelian case the exact 
calculation of the effective action is not possible, but a certain non-Abelian 
version of the Born-Infeld (NBI) action can still be envisioned in a closed form. 
This problem was extensively discussed recently ^3 E]> one reasonable 
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proposal (though perhaps still not giving the precise answer in the sense of 
superstring theory) involves a symmetrized trace ^U] construction of the action. 
Another, simpler form, uses a direct non-Abelian generalization of the U(l) BI 
action, we will call it the ordinary trace model. Investigation of static solitons 
in NBI theory ^3] has shown qualitative agreement between the results of both 
the symmetrized trace and the ordinary trace models, so here, for the reason of 
simplicity, we adopt the latter one: 



Ssbi = PJ ^x[l-^l + ^ v Ft-^{F^Fn 2 ), (1) 

with j3 being the critical BI field, in string theory (3 = X/2'Ka'. 

The simplest non-Abelian configuration for which the ordinary YM theory 
predicts chaotic behavior [H] is the following 

A = uT 1 dx + vT 2 dy, (2) 

where u and v are functions of time only, and Ti,T 2 are the gauge group 
generators. The corresponding field strength 

F = iiTxdt A dx + vT 2 dt A dy + uvT 3 dx A dy, 

contains both electric and magnetic components, but these are mutually orthog- 
onal, so the pseudoscalar invariant tr FF, generically dominant (JTJ at high field 
strength, is equal to zero. The corresponding one-dimensional lagrangian reads 

L = 2 (l - ^fl-p- 2 {ii 2 + v 2 - v 2 u 2 )^ . (3) 

The equations of motion have the form 



,2 



-uv 



2uvu(iiv + vu) 
(3 2 + u 2 v 2 



(4) 



9 2vvu(uv + vu) 



(3 2 + u 2 v 

In the limit — > oo of the ordinary YM action the second terms in the right 
hand sides of these equations vanish and one obtains the "hyperbolic billiard" (HJ 
1141 115) . known in various problems of mathematical physics, which is a two- 
dimensional mechanical system with the potential u 2 v 2 . This potential has two 
valleys along the lines u = v — 0, the particle motion being confined by 
hyperbolae uv — const. The existence of these valleys is crucial for emergency 
of chaos. 

The trajectories of the system governed by the BI action depend on the 
energy integral E and the parameter j3. Rescaling of the field variables together 
with time rescaling maps configurations with different E and /3 onto each other, 
so it will be enough to consider dependence of the dynamics only on [3 assuming 
for the energy any fixed value, e.g. E = 1. One simple (though non-rigorous) 
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method to reveal chaoticity is the analysis of the geodesic deviation equation for 
a suitably defined pseudoriemannian space whose geodesies coincide with the 
trajectories of the system ^21^7]- The lagrangian J3J gives rise to the following 
metric tensor 

ds ' 2 = ( 1 + jp\ dt2 - -jp {du2 + dy2) > (6) 

which is regular for all finite (3^0. Consider the deviation equation for two 
close geodesies 

lh ' r i^juW, (7) 



ds 

where D 2 n % is the covariant differential of their transversal separation, R l jki 
is the Riemann-Christoffel tensor, and v? is the three-speed in the space ©. 
Locally, the geodesic deviation is described by the matrix K l jkiu : >u k depending 
on points u,v,ii,v of the phase space. If at least one eigenvalue of this matrix 
is positive, the geodesies will exponentially diverge with time. Negative eigen- 
values correspond to oscillations or to slower divergence, e.g. power-law. In our 
case one of the eigenvalues of the matrix R l jkiu : >u k vanishes (due to the static 
nature of the metric), two others are the roots of a quadratic equation and read 

A 1 , 2 = i(6± v / S 2 + 4C), (8) 



where 



2uvuv — f3 2 (u 2 + v 2 ) (vu + iiv) 2 

Wuw + p 2 u 2 w ' W 

v 2 u 2 (3P 2 +u 2 v 2 ) 



p 2 u 2 w 

with W and U being positive functions 



W=l + ^ , U = l 



(3 2 ' p 2 

It is easy to see that both non-zero eigenvalues (JSJ are real, one positive and 
the other negative. Therefore, for any finite j3, there exist locally divergent 
geodesies, and it can be expected that for any /3 the motion will remain chaotic. 
But in the limit j3 — > these eigenvalues tend to zero and the analysis be- 
comes inconclusive. Qualitatively, this phenomenon of decrease of the geodesic 
divergence with decreasing (3 can be attributed to the lowering of the degree 
of chaoticity. It is also worth noting that for (3 — * the second terms in the 
equations of motion (|4I5|I look like singular friction terms. 

Another simple tool in the analysis of chaos is the calculation of the Lya- 
punov exponents defined as 

x = ^W^§ v (11) 

t^oo t \ox[0)\ 
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where Sx (t) is a solution of the linearized perturbation equation along the chosen 
trajectory. Here x stands for the (four-dimensional) phase space coordinates of 
the original system and the Cartesian metric norm is chosen. The positive 
value of x signals that the close trajectories diverge exponentially with time so 
the motion is unstable. Stable regular motion corresponds to zero Lyapunov 
exponents. 

The third standard tool is the construction of Poincare sections. Recall that 
the Poincare section is some hypersurface in the phase space. The phase space 
for the system (|13|l is four-dimensional, but the energy conservation restricts the 
motion to a three-dimensional manifold. Imposing one additional constraint will 
fix a two-dimensional surface, and we can find a set of points corresponding to 
intersection of the chosen trajectory with this surface. For the regular motion 
one typically gets some smooth curves, while in the chaotic case the intersection 
points fill finite regions on the surface. 

Both methods reveal that with decreasing (3 the degree of chaoticity is di- 
minished, though we were not able to perform calculations for very small f3 to 
decide definitely whether the chaos is stabilized or not. But in fact the ansatz 
J5J is not generic enough from the point of view of the BI dynamics, since the 
leading at high field strength term (FF) 2 in Q is zero for it. So we can hope 
to draw more definitive conclusions by exploring another ansatz for which the 
invariant (FF) is non-zero. Such an example is provided by a simple axially 
symmetric configuration of the SU (2) gauge field which is also parameterized 
by two functions of time u, v |18j : 

A = T lU dx + T 2 udy + T 3 vdz. (12) 

The field strength contains both electric and magnetic components 

F = u(T 1 dtAdx+T 2 dtAdy)+vT 3 dtAdz+u 2 T 3 dxAdy+uv(TidyAdz+T 2 dzAdx), 

and now the pseudoscalar term (FF) 2 is non-zero. Substituting this ansatz into 
the action Q we obtain the following one-dimensional lagrangian: 




u 2 (u 2 + 2v 2 ) u 2 (2iiv + vu) 2 

J 2 J 4 



L 1 =p 2 (i-Tl), n = K ji : ^ — ■ — --— ^ --(13) 



In the limit of the usual Yang-Mills theory — > oo we recover the system which 
was considered in |18j : 

L YM = \(2u 2 +i 2 -u 4 -2u 2 v 2 ). (14) 

The corresponding potential has valleys along v = 0. Like in the case of the 
hyperbolic billiard with potential v 2 u 2 , the dynamics governed by the usual 
quadratic Yang-Mills lagrangian motion exhibits chaotic character. 

However, for the NBI action with finite (3 the situation is different. One can 
use again the geodesic deviation method, the corresponding metric being: 

, 9 / u 4 + 2u 2 v 2 \ , 9 du 2 + 2dv 2 u 2 (2vdu + udv) 2 , , 

ds - I 1 + -^3^- ) M J 2 ■ (15) 
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The energy integral is given by the expression 

B + -f. (16) 

One can calculate the eigenvalues of the geodesic deviation operator as before. 
In general one finds that the larger eigenvalue is smaller than one in the case of 
the previous ansatz @. Now the regions in the configuration space appear in 
which both non-zero eigenvalues are negative, this corresponds to locally stable 
motion. With decreasing (3 the relative volume of the local stability regions 
increases, though for every j3 there still exist the regions of local instability as 
well. To obtain more definitive conclusions we performed numerical experiments 
using Lyapunov exponents and Poincare sections. The results look as follows. 
For large /3 the dynamics of the system (|13[) is qualitatively similar to that in 
the ordinary Yang-Mills theory l|14|) and exhibits typical chaotic features. The 
trajectories enter deep into the valleys along u = 0. The Poincare sections 
consist of the clouds of points filling the finite regions of the plane, while the 
Lyapunov exponents remain essentially positive. The situation changes dras- 
tically with decreasing (3. For some value of this parameter, depending on a 
particular choice of the initial conditions, one observes after a series of bifur- 
cations a clear transition to the regular motion. The points on the Poincare 
sections u — line up along the smooth curves (sections of a torus), while the 
Lyapunov exponent goes to zero. The motion becomes quasiperiodic. A typical 
feature of such regular motion is that the trajectory no more enters the poten- 
tial valleys. Some Poincare sections illustrating this chaos-order transition are 
presented on Figs. 1-3. The Lyapunov exponents corresponding to the first and 
the last picture of this series are shown on Figs. 4 and 5 respectively. 

Comparing with our first system J5J) we can guess that the substantial role in 
the chaos-order transition in the NBI theory is played by the pseudoscalar term 
in the action Since generically this term does not vanish, it is tempting to 
say that the transition observed is a generic phenomenon in the non-Abelian 
Born-Infeld theory. The effect is essentially non-perturbative in a', and we con- 
jecture that it reflects the typical smothering effect of the string non-locality on 
the usual stiff field-theoretical behavior. A related question is whether the cos- 
mological singularity, which was recently shown to be chaotic at the supergravity 
level of string models |19j , will be also regularized when higher a' corrections are 
included. Unfortunately there is no closed form effective action for the closed 
strings analogous to the Born-Infeld action for the open strings. But combin- 
ing gravity in the lowest order in a' with an exact in a' matter action we can 
probe the nature of the cosmological singularity too. This will be reported in a 
separate publication. 

The work was supported by RFBR grant 00-02-16306. 
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Figure 1: Poincare section of the sys- 
tem (JT3J) for = 0.32 



Figure 2: (3 = 0.317 
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Figure 3: /3 = 0.31 
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Figure 4: Lyapunov's exponent corre- 
sponding to Fig. ^ 
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Figure 5: Lyapunov's exponent corre- 
sponding to Fig. 
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